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We discuss entanglement and critical properties of the spin-3/2 XXZ chain in its entire gapless 
region. Employing density-matrix renormalization group calculations combined with different meth- 
ods based on level spectroscopy, correlation functions and entanglement entropies, we determine the 
sound velocity and the Luttinger parameter of the model as a function of the anisotropy parame- 
ter. Then, we focus on entanglement properties by systematically studying the behavior of Renyi 
entropies under both open and periodic boundary conditions, providing further evidence of recent 
findings about entanglement entropies of excited states in conformal field theory. 
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I. INTRODUCTION 

Since their first introduction-'^, quantum spin chains 
have been an incredibly fertile field for theoretical physi- 
cists. They are in fact interesting for various reasons: 
first of all because of their one-dimensional (ID) nature, 
that enhances the importance of quantum fluctuations 
and forbids the application of mean-field or other ordi- 
nary perturbative approaches; secondarily because of the 
integrability of some of them and the possibility of giv- 
ing a descriptions of their low-energy sector by means 
of effective quantum field theories; lastly because of the 
availability of some numerical techniques that appear to 
be particularly powerful in these cases. 

Spin- 1/2 qu antum spin chains with nearest neighbor 
interactions^i^ have been widely considered in literature 
and, among them, the XXZ spin-1/2 chain is by far the 
most studied, both analytically and numerically. Being 
integrable via Bethe Ansat^, it represents a yardstick for 
non-exact techniques. Also, the physics of its low-energy 
sector is described in the continuum li mit by a special 
class of Conformal Field Theories (CFTs)^ with confor- 
mal charge c = 1, the so called Tomonaga-Luttinger Liq- 
uids (TLLsjP^ which represent a paradigm fo r all those 
models whose excitations are of bosonic natur^SliO] 

Higher spin XXZ chains are also very interesting, be- 
cause they constitute examples of models for which, de- 
spite being non-integrable, one can provi de qu ite well 
known established field theory description^Smi More- 
over, their properties may be quantitatively determined 
with good accuracy by means of numerical simulations 
based on several efficient methods such as exact diago- 
nalization, quantum Monte Carlo, and Density-Matrix- 
Renormalization-Group (DMRG)^". Thus, they repre- 
sent a decisively indicative and efficient test in order 
to establish the validity of properties that, up to now, 
have been verified mainly in integrable systems. In ad- 
dition, the physics of these models can also be studied 
experimentally: for example, the spin-3/2 isotropic case 
is thought to model the behaviour of some kind of quasi- 
ID antiferromagnets of magnetic ions, such as CsVClaP^ 



and AgCrPsSfP, whereas various spin models may be 
now engineered in cold matter setups of trapped ional^flZl 
and have promising future implementation with ultra- 
cold atoms and molecules in strongly anisotropic optical 
lattices^^ and Rydberg atomP^T 

In recent times, much attention has been devoted 
to the connection between Quantum Phase Transi- 
tions a nd en tanglement properties in strongly correlated 
system423124j_ particular, whenever the effective theory 
describing a system is conformal, it is well known that a 
fruitful way to get physical information from numerical 
simulations, especially from DMRCt^^'^fil calculations, is 
to look at quantities known as Renyi Entanglement En- 
tropies (REs)^^ More specifically, as it will be recalled 
in the following, for TLLs the knowledge of the REs yields 
a very careful estimation not only of the central charge 
of the underlying critical theory22H2l!^ but also of the de- 
cay exponents of correlation function^2niIIIlIl^ which are 
encoded in the so called Luttinger parameter K. 

The main aim of this work is to present a complete 
investigation of critical and entanglement properties of 
the S — 3/2 XXZ model all over its gapless regime, and 
to compare the accuracy of different analysis methods 
employed to extrapolate thermodynamic quantities from 
finite-size numerical calculations based on the DMRG al- 
gorithm. In the first part, we will fully exploit the low- 
energy field theory of XXZ model and determine its rele- 
vant quantities by considering different, independent ob- 
servables; then, we will present a investigation of the en- 
tanglement entropies for bipartite intervals, considering 
both ground and excited states and systematically com- 
paring the numerical findings with CFT predictions. The 
paper is structured as follows: in Sec. [IT] we briefly re- 
view the main features of spin-S* XXZ chains from a held 
theoretical viewpoint, focusing on the TLL universality 
class emerging in the half-integer S case. We present our 
DMRG calculations and results on the 5 3/2 model 
in Sec. Ill together with a brief resume of all applied 
techniques. In Section |IV[ we perform a systematic in- 
vestigation of the Renyi entropies of both the ground and 
excited states and compare the numerical findings with 
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the predictions based on CFTs. Finally, we summarize 
the results and draw our conclusions in Sec. |Vl 



II. MODEL HAMILTONIAN 

The spin-S* anisotropic Heisenberg model, also 
known as XXZ chain, is described by the following 
HamiltoniarPI 

L 

Hxxz = Y^^S'tS^^, + SfSy^^ + A5,f (1) 
1=1 

where Si is a spin-S" operator relative to the i-th site and 
A is the anisotropy coefficient. Here, S can take positive 
half-odd-integer or integer values. In the simplest case 
S = 1/2, the model is integrable by Bethe ansatz'^, it is 
critical for |A| < 1, and, in the interval — 1 < A < 1 its 
low-energy physics is effectively described by a conformal 
field theory with central charge c = I. 

The picture becomes more puzzling as one moves away 
from the integrable S* = 1/2 case. At the isotropic point 
A = 1, chains with integer spin display a finite gap, 
whereas in the half-integer case the system is gapless and 
still described by a c = 1 CFT, as has be en prov ed in a 
series of analytical and numerical studieJ^^nHUHlSl, Jn 
such gapless regime, which persists in the finite range of 
interactions —1 < A < 1, the low-energy spectrum is 
universally described by a TLL HamiltoniarPllSl 

H=^J dx [{d^^dflK + K{d,^Y] (2) 

where d and are conjugated density and phase bosonic 
fields, while the two interaction-dependent parameters Vs 
and K are called sound velocity and Luttinger parame- 
ter respectively. Both the long-distance decay of corre- 
lation functions and spectral properties are determined 
by Vs and K: however, since no exact solution is known 
except for S — 1/2, one has to resort to unbiased nu- 
merical methods in order to quantitatively estimate the 
dependence of such parameters with respect to A. For 
the 5 = 3/2 case, that we are going to extensively study 
in the following, various numerical studies have been re- 
ported in literature, based on both exact diagonalization 
of small system^^ and simulations based on the DMRG 
algorithm^^ •^^ we will systematically refer and com- 
pare our results to the known ones in the remainder of 
the paper. 



III. NUMERICAL RESULTS: CRITICAL 
PROPERTIES 

In this section, we provide a complete study of the 
quantum critical regime — 1 < A < 1 for the S ^ 3/2 
XXZ model employing the DMRG technique, which al- 
lows to estimate both ground state and excited state 



properties with notable accuracy. Being interested in 
various physical quantities such as entropies, correlation 
functions and spectral properties, we employ simulations 
with both open (OBC) and periodic (PBC) boundary 
conditions. While the former guarantee better accuracy 
in the DMRG procedure, the latter are not affected by 
boundary effects; providing estimates in both configura- 
tion represents a good check for our final results. In order 
to accurately determine all relevant quantities of interest, 
we perform calculations with up to 512 (1156) states per 
block for OBC (PBC), together with up to 5 sweeps for 
each intermediate size during the infinite-size algorithm; 
in such a setting, typical discarded weights are of order 
10^'* (10~^) during the last iteration. 

Reliable estimates of the relevant physical quantities, 
Vs and K, may be obtained with different numerical anal- 
ysis. In the following, we will employ three alternative 
methods based on independent quantities: energy scaling 
in the low-energy spectrum, entanglement entropy and 
spin fluctuations. A detailed account on how such quan- 
tities are related to the sound velocity and the Luttinger 
parameter is given in each of the following subsections. 

A. Central charge 

As a first step in our study, we extract the central 
charge of the system from the scaling of the block Von 
Neumann Entropy (VNE), which is defined as 

Si{l,L) = -Tr APA log2 PA (3) 

Here, the system of total length L is bipartite in two 
sub-systems A,B of length l,L — I respectively, and pA 
denotes the reduced density m atrix of A with respect to 
B. In a CFT, one haPISME 

SAl,L) = ^ \og^[Lsm{7Tl/L)/TT] + s, + (4) 

where the coefficient rj — 1,2 for PBC/OBC, si is a 
model-dependent constant and 3°^'^ represents finite-size 
oscillating corrections^. In a finite system under PBC, 
the central charge can be determined by fitting the half- 
lattice entropy 

5i(L/2,L) = ^log2(L/^) + 5i (5) 

as a function of the system size L by assuming a scaling 
form of the type c{L) = cq + aoi°^ . Typical results are 
plotted in Fig. [ij a best fit of Eq. [S] for L e [28, 60] 
gives an excellent agreement with the expected value, as 
Co = 1.00 up to a 2% error over the entire parameter 
range; for negative values of A, an even better agree- 
ment is reached. Alte rnative techniques to extract c via 
finite-size scalinJ^HHl lead to comparable results. A good 
estimate of the central charge represents a reliable check 
of our numerical calculations, and is also required to con- 
sistently perform a level spectroscopy analysis without 
targeting excitations in momentum space^^ 
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FIG. 1; (Color online) Finite-size scaling of the bipar- 
tite von Neumann entropy as a function of the system size 
L £ [28, 60] for different values of the anisotropy coefficient 
A. Here, PBC are considered in order to discard oscilla- 
tory corrections. Solid lines represent fits with Eq. [5] yielding 
c = 0.995(4), 0.992(6) and 1.014(4) from top to bottom. 
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FIG. 2: (Color online) Sound velocity Vg as a function of A 
as extracted from level spectroscopy methods (black crosses); 
at A = 1 the red circle denotes the spin wave results «s = 3, 
and a previous DMRG estimatd^ Vs ~ 3.87 is indicated by 
a red star. The maximum error is of order 6 * 10~^ at the 
antiferromagnetic point. Inset: typical finite size scaling of 
the ground state energy density ega{L) for A = —0.5. 



B. Sound velocity 

Following standard level sp ectro scopy methods based 
on the energy scaling of CFTP^^ll, the sound velocity Vs 
of the single component TLL can be estimated by con- 
sidering the finite-size scaling of the ground state energy 
density, which under PBC reads as: 

eg,{L)=eo + ^ + ... (6) 

eo being the energy density in the thermodynamic limit, 
and employing the previously found values of the central 
charge. The large number of system sizes considered in 
our calculations allows to safely perform a four-parameter 
fit of the form e.gs{L) — + ai/L^ + a2/L°-^\ typical 
fitting results are shown in the inset of Fig. [2] At the 
antiferromagnetic point A = 1, our results are in good 
agreement with a previous DMRG studjEH^ where Vs was 
extracted by targeting the first excited states at finite mo- 
mentum, and in sharp disagreement with the spin wave 
resullpfil, vf^ = 3, in which quantum fluctuations are 
only approximately treated. Including logarithmic cor- 
rections according to Wess-Zumino-Novikov-WitterP^H 
theory does not lead to appreciable differences. 

The complete dependence of Vs versus A is plotted in 
Fig. [2] and suggests how, once approaching the ferromag- 
netic phase transition, the velocity of the sound excita- 
tions seems to approach zero. On the other hand, Vs 
reaches its maximum value at the Berezinskii-Koster litz- 
Thouless point^, in full analogy with the S ~ 1/2 c as 42121. 

C. Luttinger parameter 

When dealing with models whose low-energy physics 
is captured by the TLL Hamiltonian, the parameter K 



represents a fundamental quantity as it determines the 
long-distance decay of all correlation functions and thus 
which susceptibilities are the most relevant in the micro- 
scopic modeP^. We will thus employ three complemen- 
tary methods to extract such a quantity from numeri- 
cal simulations. This procedure allows to systematically 
check the validity of each method in various parameter 
regimes, and, close to the antiferromagnetic point, indi- 
cates which type of estimate is more affected by the well 
known logarithmic corrections. 

As a first way to estimate K, we make use of the pre- 
viously calculated Vg and apply the level spectroscopy 
method by targeting the first excited state with total 
magnetization Sf) = 1 under PBC. Given its en- 
ergy density e+i(L), the energy gap Ag(L) = Ljegs iL)- 
e^i{L)) scales to zero in the thermodynamic limit a d^^ l ^^ ^ 



We thus performed a least-square regression fit as a func- 
tion of l/L^ and a non- linear fit of the form Ag(L) = 
h^/L? + bi/L^^ in order to estimate the effect of higher- 
order corrections, which turns out to be negligible if 
|A| < 0.9, as can be inferred from typical finite-size scal- 
ings shown in Fig. [3] At the antiferromagnetic point, 
however, the quality of the best fit with just algebraic 
contributions turns out to be insufficient due to the pres- 
ence of strong logarithmic corrections^. We thus apply 
the same fitting procedure of Ref.[3Hl which takes into ac- 
count logarithmic corrections as ensuing from the under- 
lying SU(2) Wess-Zumino-Novikov-Witten field theoreti- 
cal structurepEH, obtaining K = 0.499 ±0.005 at the crit- 
ical point, in good agreement with previous numericaP^ 
and analytical findingj^^. A summary of the so obtained 
results all over the critical regime is presented in Fig. [8] 
and will be discussed at the end of this subsection. 
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FIG. 3: (Color online) Finite size scaling of the energy 
gap Ag as a function of for different values of A; 

solid lines are best fits (see text). The Luttinger param- 
eter estimated via Eq. [T] is, from top to bottom, K = 
2.44(0), 2.95(9), 3.69(4), 6.03(1). 



FIG. 5: (Color online) Oscillating factor of the a = 10 RE 
as a function of the system size L for OBC. The amplitude 
of the oscillations strongly decreases when approaching the 
ferromagnetic poinPl even though, an accurate estimate of 
K is still possible (with typical relative error around 10~^) by 
considering the appropriate entanglement entropy. 
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FIG. 4: (Color online) Finite size scaling of the rescaled en- 
ergy gap Ag as a function of (in logarithmic scale) for 
the same anisotropics considered in Fig. [S] solid lines are best 
fits (see text), showing how additional finite-size contributions 
are usually negligible. 



which on a finite chain of size L scales a^2ll30l35j 



(10) 



2 sin fcF— sin^ 



where kp = 7r/2 is the Fermi momentum and fa{l/L) is a 
scaling functiorP^. In order to circumvent finite size cor- 
rections to the central charge and non-universal features 
related to , one can consider the following entropy dif- 
ference: 



which, for L ^ 1, reduces tcP^. 

^4(1 + ^) 1 , cos (kpL) 



2kf 



dS^iL) 



48?7ln2fc2, L2 
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(11) 



(12) 



5„(/,i) = ^^^-^\og^[Lsm{TTl/L)/TT] + + Sr{i,L) 
brj 

(9) 

where for a — > 1 one recovers the von Neumann entropy 
^ . The influence of the Luttinger parameter on such en- 
tanglement entropies is encoded in the oscillating factor, 
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The second way of getting K is the one introduced by 
the authors in Ref. 3D] based on the finite size corrections 
of the bipartite Renyi entropies: 
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FIG. 6: (Color online) Oscillating factor of the a = 10 RE as 
a function of the system size L for PBC. Continuous lines are 
guides for the eye. 
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The finite size scaling of dSa with respect to the system 
size L allows to precisely estimates K under both OBC 
and PBC, as discussed in detail in Ref. Since it 

has been noticed that the oscillation amplitude is very 
small for Heisenberg chains with S > \/2 and a ~ 2, we 
employed REs with larger values of a in order to get a 
reliable estimate from the DMRG results. In particular, 
our estimate of K are mostly based on the a = 10 RE, as 
these data represent a good compromise between small 
oscillation amplitude and slow oscillation decajP^. 

Typical results from both PBC and OBC with sys- 
tems sizes L £ [28, 60](L € [100,180]) are presented in 
Figs. [5] and [6j We notice that even though it consider- 
ably decreases close to the ferromagnetic transition, the 
magnitude of the oscillation is still large enough to per- 
form accurate finite-size scaling. For A > 0, comparable 
results may be obtained even with smaller a's. However, 
close to the antiferromagnetic point, different types of 
corrections arise and the quality of the fitting procedure 
rapidly decreases. Thus, extracting K beyond A — 0.7 
turns out to be numerically challenging. As can be seen 
from Fig. [Sj data points obtained via dSa significantly 
deviates (of about a 5 — 15% discrepancy) from other es- 
timates in this regime. Finally, let us point out that, for 
A = —0.9, the decay exponent for PBC is so large that 
oscillations are strongly dumped for large system sizes, 
so that an accurate estimate of K is not possible; for 
OBC instead, being the exponent smaller by a factor of 
2, the estimate is still in very good agreement with all 
other methods even for A = —0.9. 

The third way we estimate K is through spin 
fluctuation^. For a bipartition as the one used to de- 
fine the von Neumann entropy, we define the quantity 

F = ((E.eA^f)') - (E.gA'g f)'- a TLL under 

PBC, spin fluctuations behave a d^" ' ^" ' "^^ !: 
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(13) 



and provide a quantitative estimation of J(I 30 | 43 | 44 ] 
employ such method by fitting the half-lattice spin fiuc- 
tuation F{L/2, L) as a function of the system size in the 
interval L g [28,60]; typical fitting results are illustrated 
in Fig.[7j At the antiferromagnetic point A = 1, logarith- 
mic corrections to correlation functions severely modify 
Eq. |13| strong oscillations, not captured by a quadratic 
TLL theory, emerge with respect to the parity of L, pre- 
venting a reliable estimate of K. 

The results obtained in this section are all shown in 
Fig. [Hj from which one can see that they agree over al- 
most the whole anisotropy parameter range. In general, 
the level spectroscopy method leads to more accurate 
estimates, as the quantities it relies on are very accu- 
rately estimated with DMRG, and additional finite size 
corrections seems to play a minor role; on the contrary, 
REs estimates, which are extremely precise for A < 0.6, 
may indeed suffer from both larger truncation errors and 
stronger finite-size dependences close to the antiferro- 
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FIG. 7: (Color online) Spin fluctuations as a function of the 
system size L under PBC. Solid lines are best fits with Eq. |13| 



magnetic point, resulting in a worst mutual agreement 
with fluctuations and LS. Finally, all results are com- 
pared with a previously stated conjecture based on ex- 
act calculations on small system size^^H -which relates the 
Luttinger parameter Ks of a spin-S' Heisenberg chain 
with 5* = 1/2 case in the A < regime: 



i^(A)5 = 25X(A)i/2. 



(14) 



Remarkably, as it can be seen from Fig.|8] this conjecture 
appears to be in a semi-quantitative agreement with the 
numerical results even well beyond its original validity 
regime. Discrepancies among the results emerge mainly 
close to the antiferromagnetic point, where logarithmic 
corrections differently affect the techniques we have pre- 
sented: in this regime, level spectroscopy turns out to 
be the most reliable method to extract the Luttinger pa- 
rameter, as fleld theoretical instruments allow to perform 
more accurate scaling hypothesis with respect to meth- 
ods based on REs and fluctuations. Nevertheless, the 
picture suggests that a deeper theoretical insight on the 
analytical properties of both REs and fluctuations close 
to critical point with logarithmic corrections may in prin- 
ciple enlarge their regime of applicability. 



IV. NUMERICAL RESULTS: ENTANGLEMENT 
ENTROPIES AND CONNECTION WITH CFTS 

A. Renyi entropies of excited states 

In recent times, an analytical formula for the REs of 
low-energy excited states in a CFT was derived and nu- 
merically verified for various quantum spin chains'*^ 
In particular, it has been predicted that the trace of the 
reduced density matrix of excited states generated 
by primary operators T of conformal weights /i, h satis- 
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FIG. 8: (Color online) Estimate of the Luttinger parameter 
as a function of the anisotropy parameter. In the legend, 
Fl, LS, RE denote spin fluctuations, level spectroscopy and 
RE results with OBC or PBC. The dashed line is the conjec- 
ture of Ref. 1381 The blue diamond (A = 0.5) and the black 
square (A = 1) indicate the LS and exact results given in 
Ref. [29| and [12] respectively. The inset shows a magnification 
of the region close to the antiferromagnetic point: here, REs 
do not provide fully reliable results due to the presence of ad- 
ditional finite-size corrections, while predictions from LS and 
spin-fluctuations are in excellent agreement. In both panels, 
numerical errors are smaller than the size of the symbols, ex- 
cept for the A = —0.9 estimate based on RE with PBC: in 
this case, the absolute numerical error due to the fitting pro- 
cedure is of order 2 due to a very strong size-dependence of 
the fits. 



fies the following relation: 



a 



(15) 



^(,) \n;ro^T(2.,7a)Tt(2.(, + 0/a);^y 



Z(l)" 



(T(0)Tt(27r/)) 



cy 



Here, a is assumed to be a positive integer, but the re- 
sult can be analytically continued all a > 1; Z{a) is 
the partition function on a torus of dimensions 2'Ka and 
2 log2[i/7r sin(7rZ/L)]; (--Ocy denotes expectation value 
on the vacuum state on a cylinder of length 27r. In par- 
ticular, if T is a vertex operator, formula ( 15 ) predicts 
that the VNE (and, up to oscillating terms, a 1 REs) of 
the excited state generated by T should be equal to the 
one of the ground state. This finding has been numer- 
ically verified in a series of exactly solvable spin mod- 
els in Ref l45|46l In order to further strengthen it, we 
considered the REs of the state generated by applying 
a vertex operator on the ground state (which belongs 



to the S' 



02 _ 



^ „j sector) thus obtaining the 
ground state in the — 1 sector. Spanning the entire 
critical region and employing PBC, we found excellent 
agreement with the CFT prediction for both von Neu- 
mann and a > 1 REs. In the former case, the entropy 
of excited and ground state coincide within numerical 
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FIG. 9: (Color online) Upper panel: von Neumann entropies 
of ground and excited states for several values of A as a func- 
tion of the block length I; here, L = 60 and PBC are consid- 
ered. Lower panel: a = 10 Renyi entropies of excited states 
a,t L — 60. Solid lines are best fits obtained using Eg. |10| with 
kp = 3l7r/60. 



uncertainty up to a constant shift of order 10~^, as can 
be seen from typical data presented in the upper panel 
of Fig. |9] for small systems up to L = 12, we further 
checked this behavior with exact diagonalization. 

We afterwards checked the relation between the oscilla- 
tion corrections and the Luttinger parameter as extracted 
in the previous section by considering REs with a > 1. 
Even in this case, the agreement with the predicted be- 
havior is remarkable, except at the ferromagnetic point 
where the quality of the fit significantly decreases. Re- 
sults of S'^[l, 60) as a function of I are plotted in Fig. [9j 
lower panel, and suggest that the amplitude of the oscil- 
lations follows a similar behavior as in the ground state 
case, namely oscillations are more pronounced when ap- 
proaching the antiferromagnetic point. 



Entanglement behavior of Heisenberg chains 
with different S 



In a recent paper, Ref.29, an interesting property con- 
cerning entanglement entropies of XXZ chains was con- 
jectured. Based on numerical DMRG simulations, it was 
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FIG. 10: A-dependence of ASi(3/2). 



shown that, under both OBC and PBC, the von Neu- 
mann entropy of XXZ spin chains with different half- 
odd-integer S > 5/2 satisfies the foUowing relation: 



(ln2)AS,iS) = S,il,L)s-S,il,L)s-i = 

(16) 

where £5 — > in the L ^ 00 limit. Such relation im- 
plies that, independently on the anisotropy A, the non- 
universal constant contribution acquires a universal form; 
interestingly, it scales as the inverse of the difference be- 
tween the one-site Hilbert space dimension of the S case 
minus the S — 1/2 one. 

We systematically investigated the behavior of 
A5'i(3/2) as a function of A by employing PBC in or- 
der to get rid of the oscillation contribution, which, as 
noticed in literatur^^l] and confirmed in the previous sec- 
tion, is extremely different for different values of S. As 
a preliminary check, we compared our value of AS'i(3/2) 
at A = 1/2 with the one reported in Refj29, finding in- 
deed very good agreement up to the different logarithmic 
basis employed here. Due to a very slight Z-dependence 
of AS'i(3/2), we mostly considered its mean value in the 
interval of block length I e [11, 29] in systems with up to 
L = 60 sites. A schematic plot of the entropy difference 
at L = 58 is presented in Fig. 10 AS'i(3/2) displays 



a notable non-monotonic A-dependence. Similar plots 
are obtained with different L. In summary, our results 



firmly confirm that Eq. 16 does not hold for S — 3/2, in 
accordance with Ref. [29l 



V. CONCLUSIONS 

In this work, we analized critical and entanglement 
properties of the spin-3/2 anisotropic XXZ chain. By 



means of DMRG simulations, we systematically calcu- 
lated the sound velocity and the Luttinger parameter 
in the entire critical region with a variety of methods 
such as level spectroscopy, entropy analysis and spin- 
fiuctuations. At the antiferromagnetic point, logarith- 
mic corrections prevent an accurate estimate of the Lut- 
tinger parameter via spin-fluctuations and entropy oscil- 
lations, whilst level spectroscopy, where such corrections 
may be systematically introduced, still provide reliable 
results. Such findings benchmark the use of level spec- 
troscopy techniques, which stem as preferable over corre- 
lation function methods based on fiuctuations and Renyi 
entropies when approaching and determining phase tran- 
sition points with logarithmic corrections. Away from 
such delicate regimes, all methods give compatible re- 
sults, in agreement with previous studies, although REs 
usually require more accurate calculations since their ab- 
solute value is very small in our case study. 

Finally, we investigated in detail the behavior of Renyi 
entropies of both ground and excited states. In the 
former case, we compared the REs of the S* = 1/2 and 
5 = 3/2 cases, proving that they are not connected by 
general relations as in the S > 3/2 case. Then, we pro- 
vided evidence that recent results on the RE of certain 
excited states in a conformal field theory are verified 
in the model of interest, presenting the first numerical 
evidence of such expected behavior in non-integrable 
models. Even though the leading corrections to REs of 
excited states and tho se of the ground state display the 
same type of behaviouil^iEll^ whether the combination of 
the two may allow for a precise determination of phys- 
ical quantities of interest still stands as an open question. 
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